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Abstract 

The protostar B335 was observed in the 1.3 mm continuum and in the H2CO 3i2 — ► 1\\ line with an angular 
resolution of about 8". The mass of the inner envelope detected by the dust continuum emission is about 0.02 M Q . 
The H 2 CO spectrum at the protostellar position shows a blue-skewed double peak profile, suggesting that the kine- 
matics of the inner envelope is dominated by infall motion. When the blueshifted and the redshifted peaks were 
imaged separately, however, there is a small east-west displacement between the maximum positions. This displace- 
ment suggests that some part of the H2CO emission might come from the outflowing gas. A combined effect of 
the infalling envelope and the outflow on the radiative transfer is discussed. This effect can make the line profile 
asymmetry severer than what is expected from infall-only models. 
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1. Introduction 

B335 is an isolated globule at a distance of 250 pc from the 
Sun (Barnard 1927 1 ; Bok et al. 1971; Bok & McCarthy 1974; 
Tomita et al. 1979; Frerking et al. 1987). In the dense core 
of B335, there is a low-luminosity (3 Lq) young stellar ob- 
ject detected in far-IR/submillimeter (Keene et al. 1983; Gee 
et al. 1985; Chandler et al. 1990). A compact source elongated 
roughly in the north-south direction was revealed by interfer- 
ometric imaging of the millimeter continuum emission from 
dust (Hirano et al. 1992; Chandler & Sargent 1993; Choi et 
al. 1999; Wilner et al. 2000). B335 is an archetypical Class 
source and considered as a protostar in the main accretion 
phase (Andre et al. 1993). It is one of the best studied pro- 
tostars. 

Early radio observations of B335 suggested that the 22 M Q 
globule contains dense molecular gas (Minn & Greenberg 
1973; Milman et al. 1975; Martin & Barrett 1978). Zhou et 
al. (1990) showed that the dense core has a steep density gra- 
dient. Star forming activities were revealed by molecular line 
observations (Frerking & Langer 1982; Frerking et al. 1987). 
CO line observations showed that there is a bipolar outflow 
oriented nearly in the plane of the sky with an opening angle 
of -50° (Hirano et al. 1988; Cabrit et al. 1988). The outflow 
axis is in the east-west direction, and the millimeter continuum 
source is located between the two outflow lobes, which sug- 
gests that the protostar at the center of the dense core is driv- 
ing the outflow (Hirano et al. 1992; Chandler & Sargent 1993; 
Choi et al. 1999). A compact radio thermal jet was detected 
near the protostellar position (Anglada et al. 1992; Reipurth 
et al. 2002). Herbig-Haro objects (HH 119) and shock-excited 
molecular hydrogen features were detected along the outflow 
axis (Reipurth et al. 1992; Hodapp 1998). 

B335 is one of the first examples of protostars showing 



1 See http://www.library.gatech.edu/barnard/ and the Dark Objects 
Catalogue therein. 



a spectroscopic signature of infall motion, i.e., blue-skewed 
double-peak line profiles. By comparing observations of sev- 
eral transitions of CS and H 2 CO with synthetic spectra from 
radiative transfer calculations of self-similar collapse models, 
parameters of the collapsing envelope, such as infall radius 
and molecular abundances, were derived (Zhou et al. 1993; 
Choi et al. 1995). More detailed observations and modelings 
using high-density tracers and near-IR extinction maps pro- 
vided a rich array of information on the protostellar collapse 
of B335 (Gregersen et al. 1997; Saito et al. 1999; Choi et al. 
1999; Wilner et al. 2000; Harvey et al. 2001; Takakuwa et al. 
2007). Interferometric observations of the millimeter contin- 
uum showed that the density gradient of the inner envelope 
has a power-law index of ~ 1 .55, which is consistent with the 
theoretical expectations (1.5) of gravitational free fall (Harvey 
et al. 2003a, 2003b). Evans et al. (2005) presented chemical 
models by comparing observations of about 30 transitions with 
model spectra generated from calculations of dust radiative 
transfer, molecular line radiative transfer, and chemical evo- 
lutions. These observations and modelings collectively showed 
that B335 is probably the best object for the study of proto- 
stellar collapse and related phenomena, and it is important to 
collect more data in order to test our understanding of star for- 
mation at the earliest evolutionary stage. 

This Letter presents observations of the B335 region in the 
A = 1.3 mm continuum and the H2CO line using the Berkeley- 
Illinois-Maryland Association (BIMA) array. In § 2 we de- 
scribe our observations. In § 3 we report the results of the 
BIMA imaging and discuss the structure of the protostellar 
core. 

2. Observations 

The A = 1.3 mm continuum and the H2CO 3i2 — > 1\\ line 
(225.697772 GHz) were observed using the BIMA array in the 
D-array configuration on 1999 September 4. The phase track- 
ing center was a 200 o = 19 h 37 m 00?99 and <5 2 ooo = 07°34'10."8. 
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A narrow spectral window was used for the H2CO line with a 
spectral resolution of 0.098 MHz giving a velocity resolution of 
0.13 km s -1 . For the continuum, six wide spectral windows at 
each sideband were used, resulting in a single-sideband band- 
width of 600 MHz. These wide windows do not contain any 
strong line. The phase was determined by observing a nearby 
quasar 2013+370 (NVSS J201528+371059). The flux calibra- 
tion was done by observing Uranus. Maps were made using a 
CLEAN algorithm. Natural weighting was used for all the im- 
ages presented in this Letter. 

3. Results and Discussion 

3.1. Millimeter Continuum 

The 1.3 mm continuum emission was detected in both the 
lower and the upper sidebands. Since combining the data of 
the two sidebands did not improve the signal-to-noise ratio, 
we present the map made using the lower sideband data only 
(Fig. 1). The peak flux density is 97 ± 19 mJy beam" 1 , and 
the total flux is 180 ± 60 mJy. (The peak flux of the upper- 
sideband map is 70 ± 20 mJy beam -1 .) The 1.3 mm peak po- 
sition is consistent with the millimeter/centimeter position of 
previously published maps (Wilner et al. 2000; Reipurth et al. 
2002). 

Figure 2 shows the spectral energy distribution of B335 in 
the millimeter-centimeter range. From the three data points 
in the millimeter band, the best-fit power-law spectrum has a 
spectral index of a « 2.3, which suggests that the millimeter 
emission mostly comes from dust. For comparison, the spectral 
index in the submillimeter band within a 40" aperture is 2.9 ± 
0.6 (Shirley et al. 2000). 

To derive the mass from the dust continuum flux, the dust 
emissivity given by Beckwith & Sargent (1991) is assumed, 
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Fig. 1. Map of the A = 1.3 mm continuum toward the B335 region. 
The contour levels are 40 and 80 mJy beam -1 . Dashed contours are 
for negative levels. The rms noise is 19 mJy beam -1 . Shown at the 
bottom left corner is the synthesized beam: FWHM = 10."6 X 6."7 and 
P. A. = 3 1 ° . The straight line at the bottom corresponds to 2000 AU at 
a distance of 250 pc. Plus: Position of the 1.2 mm continuum source 
(Wilner et al. 2000). Filled circle: Position of the 3.6 cm continuum 
source (Reipurth et al. 2002). 
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Fig. 2. Spectral energy distribution of B335 from interferometric ob- 
servations. Flux densities are from Reipurth et al. (2002), Choi et al. 
(1999), Chandler & Sargent (1993), and this work. The flux density at 
A = 3.4 mm from the data presented by Choi et al. (1999) is 21 ± 4 
mJy. Dotted straight line: Best-fit power-law spectrum in the millime- 
ter range. The spectral index is a = 2.3 ± 0.3. 

/ \ 

K, y = 0.l(^j cm 2 g-\ (1) 

where v is the frequency, i/q = 1200 GHz, and (3 is the opacity 
index (f3 m a — 2, in the millimeter range). If the emission is 
optically thin, the mass can be estimated by 

F D 2 

M=^— (2) 

where F v is the flux density, D is the distance to the source, 
B v is the Planck function, and Td is the dust temperature. 
Dust radiative transfer modelings show that the dust temper- 
ature within 0.006 pc (5") from the center ranges from 20 K 
to -40 K (Evans et al. 2005). Assuming D = 250 pc and T d 
= 30 ± 10 K, a rough estimate of the mass from the interfer- 
ometric observations is M = O.OH^q qq^ Mq, which includes 
the inner protostellar envelope and the disk. The mass estimate 
is very sensitive to the opacity index, and Chandler & Sargent 
(1993) derived 0.21 M Q by assuming (3= 1.2. For comparison, 
the 1.3 mm flux density within a 40" beam is 0.57 Jy, and the 
mass estimate within a 120" beam is 1.2 Mq (Shirley et al. 
2000). 

3.2. The HiCO Line 

Figure 3 shows the map of the H2CO line integrated over 
the full width of the line. The strongest emission peak co- 
incides with the central protostellar position. There is a sec- 
ondary emission structure near the northwestern corner of the 
map, which is probably related with the outflow. 

Figure 4 shows the H2CO spectrum at the protostellar po- 
sition. Similarly to most of the optically thick molecular lines 
seen toward B335, the H2CO line also shows a blue-skewed 
double peak line profile that may be interpreted as a sign of in- 
fall motion. However, when the blueshifted and the redshifted 
peaks were imaged separately, there is a small east-west dis- 
placement between the maximum positions (Fig. 5). The spec- 
tra at the maximum positions are shown in Figure 6. Since the 
bipolar outflow is also in the east-west direction (Hirano et al. 
1992), the displacement may be related with the outflow. That 
is, some part of the H2CO emission might come from the out- 
flowing gas. This contribution from the outflow may provide 
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Fig. 3. Map of the H 2 CO line toward the B335 region. The H 2 CO 
line intensity was averaged over the velocity interval of Vlsr. = (7.6, 
9.0) km s -1 . The contour levels are 1, 2, and 3 times 0.7 Jy beam - 1 . 
The rms noise is 0.26 Jy beam -1 . Shown at the bottom left corner is 
the synthesized beam: FWHM = W'.'6 X 6'.'1 and P.A. = 32°. Plus: 
Position of the 1.2 mm continuum source (Wilner et al. 2000). 
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Fig. 4. Spectrum of the H2CO line at the B335 protostellar position. 
Vertical dotted line: Systemic velocity of the B335 cloud core, based 
on the average of optically thin lines (Vlsr = 8.30 ± 0.05 km s ; 
Evans et al. 2005). 

an explanation on the particularly poor fit of the H2CO 3 12 — > 
2n line in the modeling of line profiles (Zhou et al. 1993; Choi 
et al. 1995). 

Wilner et al. (2000) arrived at a similar conclusion from 
interferometric observations of the CS J = 5 — > 4 line. The 
high-resolution CS maps are dominated by clumps aligned in 
the east-west direction, which was interpreted as a slow, dense 
component of the bipolar outflow. However, there is an interest- 
ing difference between the H2CO and the CS maps. The strong 
blueshifted peak of the H2CO emission is located ~3" east of 
the protostar (Fig. 5), while the strong blueshifted CS emission 
peak is located ~5" southwest (Wilner et al. 2000). If the blue 
asymmetry of the line profile was caused entirely by a bright 
clump in the outflowing gas, the strongest blueshifted emission 
should have come from the same part of the outflow. The differ- 
ence between the H2CO and the CS maps suggests that the out- 
flow alone is not necessarily the dominating cause of the line 
profile asymmetry. Since most of the optically thick molecu- 



Fig. 5. Maps of the blueshifted (dashed contours) and the redshifted 
(solid contours ) peaks of the H2CO line. The H2CO line intensity was 
averaged over the velocity intervals of Vlsr = (7.7, 8.2) and (8.4, 8.9) 
km s , respectively. The contour levels are 1, 2, and 3 times 1.3 Jy 
beam -1 . The rms noise is 0.44 Jy beam -1 . 
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Fig. 6. Spectra of the H2CO line at the maximum positions of the 
blueshifted (top) and the redshifted (bottom) emission. The position 
offset in arcseconds relative to the map center is written at the top left 
comer of each panel (see Fig. 5). 

lar lines of B335 show a blue-skewed line profile (Evans et 
al. 2005), the line asymmetry observed with single-dish tele- 
scopes is probably caused by a systematic effect, most likely 
the collapse, rather than the clumps in the outflow. Especially, 
the blue-skewed profiles of off-center spectra at positions in 
the north or south of the protostar (i.e., at positions in the di- 
rection perpendicular to the outflow) provide a strong evidence 
for infall motion (Zhou 1995; Choi et al. 1999). 

While the difference between the CS and the H2CO maps 
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may be caused by the chemical differentiation, the H2CO and 
the CS maps have two features in common. First, the overall 
shape of the emission structure is elongated in the east-west 
direction, parallel to the outflow axis. Second, the blueshifted 
emission is either stronger than or comparable to the redshifted 
one in both the eastern and the western sides of the protostar. 
One possible explanation of these features may be the com- 
bined effect of the infalling envelope and the outflow on the 
radiative transfer (Fig. 7). If the outflow axis is nearly in the 
plane of the sky, and if the outflowing gas and the infalling gas 
along a single line of sight have comparable (projected) veloc- 
ities, the emission from the redshifted outflowing gas can be 
absorbed by the redshifted infalling gas in the front hemisphere 
of the envelope. In contrast, the emission from the blueshifted 
envelope may not suffer from absorption if the blueshifted out- 
flowing gas is optically thin and/or more highly excited. This 
effect can make the line profile asymmetry severer than what is 
expected from infall-only models. This model may not apply to 
every species/transitions because the infalling and the outflow- 
ing gas have quite different physical and chemical conditions. 
For example, the enhancement of gas-phase H2CO and CS 
is expected in shocked outflowing gas while some molecules 
(such as DCO + ) are abundant only in the quiescent gas (for 
example, Bachiller & Perez Gutierrez 1997; Blake et al. 1995). 
Such an outflow-envelope coupling in radiative transfer may 
work for high-density tracer lines that are optically thin in the 
outflowing gas and moderately optically thick in the infalling 
envelope. For a better understanding and comparisons with ob- 
servations, this model needs to be tested in detail using multi- 
dimensional radiative transfer codes. 

We thank KASI ISMJC members for helpful discussions. 
This work was partially supported by the LRG Program of 
KASI. 
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